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Q-i! We present an analytical calculation of the leading three-loop radiative 

O 1 correction to the 5-parameter in the Standard Model in the large Higgs mass 

limit. Numerically, = 1.1105 x g 4 /(10247r 3 ) x mjj/Myy. When combined 
with the corresponding three-loop correction to the p-parameter, this leads to 



shifts of A^ 3 ) sin 2 O^r^ = 4.6 x 10 9 x mj^/M^ in the effective weak mixing 

^ ! angle and A^Mw = —6.3 x 10~ 4 MeV x mjj/Myy in the W boson mass. For 

both of these observables, the sign of the three-loop correction is equal to that 
of the one-loop correction. 
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1 Introduction 



Radiative corrections to four-fermion processes are sensitive to the presence of heavy 
particles in the standard model. This allows one to draw conclusions about the 
mass of the Higgs boson from fits of precise electroweak measurements to theoreti- 
cal predictions of the standard model |TJ. For a light Higgs boson, the Higgs mass 
dependence of the theoretical predictions is mainly due to one-loop radiative cor- 
rections to the gauge boson propagators. These one-loop corrections depend on tuh 
logarithmically [2 j. However, because the Higgs self-interaction A is proportional 
to m 2 H , there are higher order radiative corrections which grow like powers of m# 
in the limit m# — > oo, eventually overcoming their relative suppression by powers 
of a. Such higher order corrections could become important if the Higgs boson is 
very heavy. At the two- loop level, the leading corrections are proportional to m#, 
but the numerical coefficient of these terms turns out to be very small [13 IH E] , and 
therefore, they are not important for tuh less than a few TeV. However, it has been 
suggested that the smallness of the two-loop corrections may be somewhat acciden- 
tal jSj. If this is true, then one may expect larger corrections to appear at the three 
loop level. 

An explicit calculation of the leading three-loop correction to the electroweak p- 
parameter in the large Higgs mass limit [7] has shown that, for this observable, the 
numerical coefficient of the three-loop correction is indeed larger than the two-loop 
coefficient, and that the leading three-loop correction Ap^ is already equal in mag- 
nitude to the two-loop correction Ap^ at a Higgs mass of approximately 480 GeV. 
The purpose of the present paper is to extend our previous investigation to other 
electroweak observables, in particular the sine of the effective leptonic weak mixing 

angle sin 2 O^r^ , defined in terms of the couplings of the Z-boson to leptons, and 
the mass of the W boson. At the two-loop level, the complete electroweak fermionic 

corrections to sin 2 are known |SJ , and for the W mass, both the fermionic and 
the bosonic corrections have been calculated p. Here, we are concerned with the 
leading three-loop bosonic corrections, which grow like m# in the large Higgs mass 
limit. 

We will present our results in terms of the S, T and U parameters introduced by 
Peskin and Takeuchi ^U] to describe the effects of heavy particles that enter only 
through corrections to the gauge boson propagators (oblique corrections). They 
are defined in terms of the transverse gauge boson self-energies at zero momentum 
Yj* = S^(p 2 )| p 2 =0 and their first derivative = (p 2 )| P 2 =o as 
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where cw = cos 6\y, sw — sin 6^. The leading corrections to four-fermion processes 
in the heavy Higgs mass limit can be treated in this framework provided that there 
are no further contributions from vertex and box diagrams. In the renormalization 
scheme we are using, where m# and Mw are defined by subtractions at p 2 = —m 2 H 
and p 2 = 0, respectively (which is equivalent to the on-shell scheme if one is only 
interested in the leading term in the large Higgs mass limit), this condition is satis- 
fied EE]. 

2 Calculation 

All our diagrams are generated by the program QGRAF ^2]- The rest of the 
calculation is done mainly in FORM ^3]. First of all, we need the bare gauge boson 
self-energies, which we decompose into their transverse and longitudinal components 

as 

SJ(P) = " P ~f) ^(P 2 ) + V ~f- ^W) , (4) 

where X = AA, AZ, ZZ,WW . The scalar functions S^(p 2 ) and S^(p 2 ) are ex- 
tracted from by means of the projectors Pt = (d\w — PuPu/p 2 ) and 
Pl — PfiPv/p 2 respectively, where d = 4 — e is the space-time dimension, and ex- 
panded in a Taylor series up to order p 2 . The coefficients are vacuum integrals, 
which, in general, depend on three different non-zero masses: m#, Mw and Mz- 
We use the method of expansion by regions ^3] to expand them in powers of van for 
tuh 3> Mw,Mz, keeping all terms proportional to or higher in the three- loop 
gauge boson self-energies. As discussed in more detail in ref. [7j, the expansion by 
regions produces: 

(a) three-loop single-scale vacuum integrals, where the mass scale is mn (from the 
region where all loop momenta are large); 

(b) products of one- (two-) loop vacuum integrals depending on mn and two- (one-) 
loop vacuum integrals depending on the small scales Mw and Mz (from the regions 
where some of the loop momenta are small and others are large); 

(c) and finally three-loop vacuum integrals depending on Mw and Mz (from the 
region where all loop momenta are small). 

Using integration-by-parts identities, we reduce all of these vacuum integrals to 
master integrals. For cases (a) and (c), we use the automatic integral reduction 
package AIR [TS]. The integrals that occur in case (a) have Af + < 5 and M- < 2, 
where Af + is equal to the sum of the positive indices, minus the number of positive 
indices, and M- denotes the absolute value of the sum of the negative indices. 
The three-loop single-scale master integrals appearing in case (a) can be found as 
expansions in e in ref. |16| ITTj . 

In case (c), we have to reduce two-scale integrals with A/+ < 4 and Af- < 2. 
However, once we sum over all diagrams, in all the three-loop gauge boson self- 
energies, the non-factorizable three-loop master integrals occurring in case (c) cancel, 
so that explicit formulae for such master integrals are not needed. 
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The longitudinal parts of the gauge boson self-energies are related to the self- 
energies of the scalars and the mixings between the scalars and the gauge bosons 
by a set of Ward identities, which are shown explicitly in appendix A. We have 
verified that these Ward identities are satisfied by the full (ie., including tadpole 
contributions), unrenormalized self-energies (up to order p 2 ). 

3 Renormalization 

We work with the Lagrangian, C = C inv + Cf ix + Cpp, with the invariant part 
C inv = ~w; u W a ^ - - ( J D M $) t - iA($ t $) 2 (5) 



with the Higgs doublet 



and the gauge fixing term 

C flx = -C+C-- l -{C z f- l -{C A )\ (7) 
with 

C ± = -d^W** + M W (j) ± (8) 
C z = -d^Z" + M z <j>° (9) 
C A = -dpA". (10) 
The parameters v, A and fi are given by: 

. = V2^, (11) 
9 



x - * 2 Sr- (12) 



2 
H_ 

2 

W 



ix = (3- l -m 2 H . (13) 

Giving (3 a non-zero value produces a term in the Lagrangian that is linear in the 
Higgs field H . We adjust this parameter order by order to make the renormalized 
Higgs tadpole vanish. We introduce two further renormalization constants, Zh and 
Z mH , by making the following substitutions in C inv (but not in Cf ix or Cfp)' 

Zm H 

m H (14) 
M w -> Z H M W (15) 
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H — > Zh H (16) 
0± ^ (17) 

0° - W. (18) 

The renormalization constants are fixed by imposing the following conditions on the 
renormalized Higgs tadpole T H,ren and on the renormalized <j) and H self-energies: 



p(l),-ff,ren 


= o«) 


(19) 


-p(2),H,ren 


= c«) 


(20) 


L | p2=0 


= o«) 


(21) 


L | p 2 =0 


= o«) 


(22) 


p„v(l)^™ 

rve zj |p2 +m 2 f =0 


= oW) 


(23) 


■ n ' e ^ \p 2 +m 2 H =0 


= 0«) 


(24) 



These renormalizations remove all the terms of order m 2 H and of order mjj from 
the one- and two-loop gauge boson self-energies, the self-energies, and the mixings 
between 0's and gauge bosons. This ensures that no two- or three-loop vertex or box 
graphs containing such self-energies as subgraphs can give corrections that grow like 
m 2 H or run in the large Higgs mass limit, and allows us to restrict our attention to 
the corrections that come from the gauge boson self-energies. The renormalization 
procedure used here is slightly different from the one used in ref. where the 
parameter f3 was not used and therefore, tadpole contributions had to be taken into 
account explicitly. In particular, the constant Z niH used here is different from Z m of 
ref. |7j . However, the two procedures are equivalent to each other. Explicit formulae 
for the renormalization constants are listed in appendix B. 

We have verified that the renormalized longitudinal photon self-energy and photon- 
Z mixing are zero, as they should. 



4 Results and conclusion 

Here, we collect the leading one-loop [TBI EH] and two-loop [H |H] contributions and 
the new three-loop contribution to the ^-parameter in the heavy Higgs limit: 

S {1) = T^logf^ , (25) 



4tt V 16tt 2 J A!*, V 72 8 54 
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4tt V 16tt 2 / 



10 3456 " v ' 24 1* 
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The constants appearing in these expressions are ^01 \H\ : 



E/ 3 ,i = ^C( 4 ) + ^C(2)log 2 2-^log 4 2-Li 4 ^) = -0.117S75999G5C2S) 



-0.03901272636, (29) 



Vs.1 
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m 3 n 



1.0149416064, (30) 



log 2 



2 sin — 
2 



-2.1447672126. (31) 



The T-parameter is related to the p-parameter by Ap = aT, so that, from refs. [23 

013, 

^ = -rAr^(^) > (32) 



16^ b \M^ 



T(1> = ^\-&^M^ m) ' (33) 

Ti3 > - M^) 2 M ( - U282h m 

In contrast to S, the parameter U can only be different from zero if custodial sym- 
metry is broken, and as a result, it should be suppressed compared to S 5J. In the 
approximation where we keep only logarithmic terms in m# at the one-loop level, 
quadratic terms at the two-loop level and quartic terms at the three-loop level, we 
find that U vanishes. This provides a very useful check on the calculation. 

Knowing S, T, and U, one can easily obtain the corresponding corrections to 
various electroweak precision observables fUJ 122]- We illustrate this by two 
examples. The effective weak mixing angle is shifted relative to its tree level value, 
expressed in terms of a, Gp and M z , by 



sin 2^1ept_ A 2 ^lept 1_ I Tra 
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en 
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en 
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3.8 x 10" 4 


-6.7 x 10~ 8 


7.4 x 10" 8 


5 


8.9 x 10" 4 


-4.2 x lO" 7 


2.9 x lO" 6 


10 


1.3 x 10~ 3 


-1.7 x lO" 6 


4.6 x 10~ 5 


15 


1.5 x 10~ 3 


-3.8 x 10~ 6 


2.3 x lO" 4 


20 


1.6 x 10~ 3 


-6.7 x 10~ 6 


7.4 x 10" 4 


25 


1.8 x 10- 3 


-1.1 x 10~ 5 


1.8 x 10~ 3 



Table 1: Corrections to sin 2 as a function of m^/Mjy 



m H /M w 








2 


-0.055 


0.000041 


-0.000010 


5 


-0.13 


0.00025 


-0.00039 


10 


-0.18 


0.0010 


-0.0063 


15 


-0.21 


0.0023 


-0.032 


20 


-0.24 


0.0041 


-0.10 


25 


-0.26 


0.0064 


-0.25 



Table 2: Corrections to My/ in GeV as a function of ttih/Mw- 



with 



A sin 2 6 



lept 
eff 



a 



r 2 - s 2 



Similarly, one finds a shift in the W-mass, 



(36) 



M w = AM W + M z 



\ 



+ 



with 



A Mi 



w 



w 



V2G F Mj 



S + c 2 w T + 



C W S W jj 



Numerical values of these shifts, calculated using 

n e 2 Air a 



9 = ^ 



(37) 



(3* 



(39) 



for the weak coupling constant, with a = 1/137 and s 2 ^ = 0.23, are shown in Tables^ 
and El and in Figures ^ and El The two- loop corrections are extremely small. This 
is partly due to a cancellation between and in eqs. (J3TI|) and (|3*5|). At the 
three-loop level, there is no such cancellation between and , since they have 
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opposite signs. For both observables, the three-loop correction becomes equal to 
the one-loop correction at m# ~ 2 TeV. Since the three-loop term has the same 
sign as the one-loop term in both cases, it is clearly not possible for the radiative 
corrections due to a heavy Higgs boson to mimic those of a light Higgs boson. 
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A Ward identities 

The Ward identities for the self-energies can be derived by expanding the following 
identities for the dressed propagators G: 



1 


= ffG A ^{v) 


(40) 





= P»P u G* u z (p)+t P »M z Gf A (p) 


(41) 


1 


= P^G^ip) + 2ip»M z Gf z (p) + M 2 z G^\p) 


(42) 


1 


= P»P v G w ^{p) + 2ip»M w Gf v (p) + M 2 w G^{p) 


(43) 



Expanding the propagators G in eqs. (}4T)|) - (j4*3*j) in terms of the self-energies S, one 
obtains a set of Ward identities. For the mixings between Goldstone and gauge 
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Figure 2: Shifts in M w as a function of m H /M w . 

bosons, we write 

where X = <j> Q A, (j)°Z, <j>W . 

Up to three loop order, we find the following relations. 

= E^ AA 

= ^ AZ + iM z E^° A 

= p 2 ^' zz + 2iM zP 2 ^° z + MlE«^° 

= p 2 Y$' ww + 2iM wP W>>* w + M^E«^ 

= e«^ 1 -(eW.^) 2 

o = ^'^-^(s^ ^ 

= p 2 Ei 2) ' w - p 2 (E«^) 2 + 2 zM w (E^^E^^ + ^E^) 
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AW 
A« + A( 2 ) 
A« + A( 2 ) + A (3) 



+ M^^(E«^) 2 + E( 2 )^ (52) 
= E P).^ _ i E (i),^V° ^(i),fij 2 _ 2 s (i),«°A E (2),^4 (53) 

= S (3),AZ _ j_ E (l)^V S (l)^°A E (l)^°Z _ S (1),^A S (2),0°Z _ S (1),^Z E (2),^A 
+ iM z ( 1 f 2 S (1),0°A + ^ S (1),0V O S (2),0 O A + ^ S (2)>V° S (1),0 ^ 

+ £( 3 ^° A ) (54) 
= p 2 t!1 ] ' zz - (e (1) ^° z ) 2 - 2 p 2 S (i)^ o z s (2),0°z 

+ 2iM z ^J_^£(iWVj 2 s (i),<^ + s (i)^ ( V° s (2)^°z + S (2),^0 O S (1),0 ^ 

+ p 2 E (3),^ 

+ Mf ( J^E^ ) 3 + ^ S (i)^V° s (2)^V + E (3),^ (55) 
= p 2 y!1 ] ' ww - yj(i)><«>(yj(i)><^) 2 - 2p 2 S (1) '^S (2),0W 

+ 2 iM w f (E (1) '^) 2 S (1) '^ W + sW^C 2 )^ + s( 2 )'^S (1) '^ 



+ p 2 E (3),^^ 

B Renormalization constants 

In d = 4 — £ dimensions, the renormalization constants Zh and (3 are given by 
Z H = 1 - <S£> - 5g } and /3 = m 2 , + /?( 2 )) , with 

= uki'^W^T) (57) 



(2) 



1 <7 



2 



i(27r) d M2, 



A«;1) 2 (4-t7^t+ ' 



128 16 (e - 4) 8 (e - 4)< 



3s 9 



+ m 2 H I 2 (m 2 H , m 2 H , m 2 H ; 1, 1, 1) ( ^ + ^ + 
+ m 2 H I 2 (m 2 H: 0,0; 1,1,1) 



64 64 8 (e - 4) 
£ 11 3 



64 64 4 (e - 4) 



+ — 7+0 + 



i(2vr) d M2 



w 



,6™ 

4 8 (£-4),' mff 



/3« 



(58) 
(59) 



1 4 2 



A«; 1)M-^- 



32 8 (5 - 4) 



15 9e 



+ m 2 H I 2 (m 2 H , m 2 H , m 2 H ; 1, 1, 1) ( £f - ^ 

+ / 2 K,o.o;M,i)K 



64 64 
9-3e 



2 , ... M2 (2c 2 w + l)(e-2)( £ -9) 



64 



w 



32 c 2 ^ 



+ h{m 2 H - l)h{M 2 w - 1) ['?£ - A) + l^Mf; 1) (?£ - 1 



+ 



1 £ 2 
^(27r) d 



J iK; i) (| - y) + i) + ^i(Mf; i)} e 



(60) 



The one- and two-loop scalar integrals used in these expressions are defined by 
I 2 (m\, ml, ml ; n u n 2 , n 3 ) = 

J d% d% P(h ; m 1 ) ni P(k 2 ; m 2 ) ri2 P(A; 1 + k 2 ; m 3 )" 3 (61) 



Ji(m 2 ;m) 



/ 



d d fci ; mi) 



rai 



with 



P(/c ; m) = 



1 



(62) 



(63) 



k 2 + m 2 

Note that in (3^ and /9^ 2 \ it is necessary to keep not only the leading power of m H , 
but also the next-to-leading power. 
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The Higgs mass renormalization constant Z m[I = 1 — 5m H — Sm H 

V 32 16 16 8 2/ J ' 

A(2) ( ? \*( m *Y t j*(<i e \ m sj 27 V 27 /5 363 



(64) 



1575 27 



256 64 



64 & 128 32 512 32 32 sv ; J 



(65) 
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